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THE  EFFECTS  OF  NEUTRAL  GAS  RELEASE  ON  VEHICLE  CHARGING: 

EXPERIMENT  AND  THEORY 


I.  Introduction 

As  is  well-known,  spacecraft  in  LEO  can  acquire  large  negative  potentials  relative  to  the 
local  plasma,  and  these  potentials  can  damage  sensitive  spacecraft  electronics  even  in  the  absence 
of  active  experimentation  as,  for  example,  in  the  case  of  solar  arrays  [  McPherson  et  al.,1976; 
Hastings  et  al.,  1992 J.  Such  large  negative  voltages  can  be  acquired  by  movement  through 
natural  environments  such  as  the  auroral  regions  [Gussenhoven  et  al.,  1985;  Mullen  et  al,  1986], 
or  because  of  active  space  experiments  which  emit  a  charged  particle  beam  when  insufficient 
neutralizing  return  current  is  collected  [Borovsky,  1988].  For  example,  at  negative  charging 
levels  the  low  mobility  of  positive  ions  impedes  their  ability  to  act  as  neutralizing  return  current. 
Another  source  of  payload  charging  is  ion  beam  emission  [Kaufmann  et  al.,  1989;  Roth  et  al., 
1983] ,  but  typically  the  currents  are  small  and  charging  has  not  been  observed  above  a  few 
hundred  volts  [Olsen  et  al.,  1990].  However,  electron  beam  sources  employed  in  many  space 
experiments  are  capable  of  relatively  large  currents  and  they  can  also  cause  high-voltage 
charging;  both  positive  and  negative  charging  is  predicted  to  occur  because  of  ringing 
[Borovsky,  1988;  Winglee,  1991].  Also,  satellites  in  eclipse  in  the  auroral  region,  for  example, 
and  geosynchronous  spacecraft  can  often  reach  high  negative  voltage  levels.  These  voltages  can 
reach  hundreds  of  volts  and  have  been  observed  as  high  as  several  kilovolts.  In  addition, 
planned  charging  and  discharging  events  typically  produce  a  rapid  change  in  the  payload 
potential  with  a  slow  decay  to  equilibrium  so  that  large  potentials  can  occur,  for  example,  on  the 
LEO  spacecraft  solar  arrays  mentioned  above  [Hastings  etal,  1992].  Conventionally,  most 
solar  array  systems  for  U.S.  space  vehicles  have  a  bias  of  28  Volts.  Future  solar  arrays,  however, 
are  being  designed  for  much  higher  voltages  to  meet  higher  power  demands  at  low  currents.  The 
space  station,  for  example,  has  160  Volt  arrays.  Resulting  electrical  discharges  in  solar  arrays 
have  led  to  performance  losses  in  several  recent  satellite  launches. 

Because  of  the  above  considerations  and  the  active  nature  of  desired  experimentation, 
much  of  the  focus  of  the  Space  Power  Experiments  Aboard  Rockets  (SPEAR)  m  program  was 
the  study  of  a  variety  of  spacecraft  grounding  techniques.  This  was  part  of  the  overall  objective 
of  investigating  the  interaction  of  high  voltage  power  systems  in  the  low  earth  orbit  (LEO) 
environment.  The  SPEAR  program  was  begun  by  the  Ballistic  Missile  Defense  Organization 
(BMDO)  in  the  late  1980's  and  continued  with  the  successful  launches  of  two  sounding  rockets, 
in  addition  to  a  series  of  ground  tests  in  large  laboratory  environments  [Berg  et  al.,  1995; 

Gilchrist  et  al.,  1993;  Greaves  etal,  1990;  Alport  et  al.,  1990].  Several  of  the  results  from  these 
tests  showed  promising  possibilities  for  mitigation  of  the  charging  through  neutral  gas  release 
[Berg  et  al.,  1995].  An  attractive  feature  of  this  method  lies  in  the  fact  that  most  spacecraft  have 
relatively  large  amounts  of  neutral  gas  storage,  often  for  the  purpose  of  attitude  control 
adjustments.  Therefore,  if  this  reservoir  can  be  used  in  mitigation  techniques  it  might  provide  a 
means  of  ameliorating  the  charging  without  the  incorporation  of  specially  designed  systems  (such 
as  plasma  thrusters).  Recent  results  [Mandel  et  al.,  1998;  Machuzak  et  al.,  1996;  Berg  et  al., 

1995;  Gilchrist;  1993;  1990]  suggest  that  this  approach  may  hold  promise  in  typical  operational 
scenarios,  both  in  rocket  flights  and  Space  Shuttle  configurations. 
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This  paper  describes  laboratory  experiments  designed  to  investigate  the  effects  of  Neutral 
Gas  Release  (NGR)  on  space  vehicle  charging.  The  effort  is  an  outgrowth  of  the  SPEAR 
program.  The  experimentation  includes  use  of  a  simulated  “payload”  cylinder  which  is  situated 
in  a  large  vacuum  chamber  at  the  Naval  Research  Laboratory  (NRL)  and  which  incorporates  a 
neutral  gas  release  valve.  The  cylinder  and  valve  are  charged  to  kilovolt  levels  (usually 
negative)  with  respect  to  the  chamber  walls,  both  in  vacuum  and  in  plasma  environments,  and  the 
neutral  gas  is  released,  resulting  in  discharge,  or  breakdown,  to  the  local  plasma.  The  charging 
electronics  package,  which  was  built  at  NRL,  was  used  in  the  SPEAR  laboratory  testing  and  was 
afterward  refitted  for  use  in  the  Space  Physics  Simulation  Chamber  (SPSC).  The  experiment 
consists  of  varying  the  mass  of  the  release  gas,  composition,  flow  rate  and  mach  number  of  the 
releases  and  determining  the  effects  on  final  voltage  and  voltage  discharge  initiation  levels. 

The  theoretical  treatment  of  the  discharge  is  treated  through  an  analysis  which  employs 
the  Townsend  discharge  criterion.  The  fundamental  assumption  in  this  model  is  that  when  a 
negatively  charged  object  in  neutral  gas  discharges  (or  when  the  discharge  takes  place  from  the 
cathode  to  the  anode)  the  number  of  electrons  grows  exponentially  in  an  avalanche  fashion  along 
the  discharge  path.  The  model  then  is  highly  neutral  gas  dependent  through  the  Townsend 
ionization  coefficient  which  is  defined  as  the  number  of  ionization  events  occurring  per  unit 
length  in  the  electric  field  direction.  In  addition,  the  nozzle  shape  and  design  have  a  fundamental 
effect  on  the  neutral  gas  distribution  and  therefore  on  the  discharge  characteristics.  We  treat 
these  areas  in  a  simple  model  which,  through  flux  conservation  and  the  assumption  of  lsentropic 
expansion  upon  release,  appears  to  successfully  model  the  physics  of  the  release  and  its  relation 

to  the  subsequent  discharge. 

II.  Experimental  Description 

The  experimentation  was  performed  in  the  Naval  Research  Laboratory’s  Space  Physics 
Simulation  Chamber  (SPSC)  using  an  existing  diagnostic  complement  developed  for  this  and 
other  programs.  One  particular  advantage  of  this  experimental  environment  is  the  large  size  of 
the  facility  (1.8  m  x  5  m  vacuum  chamber)  in  addition  to  the  ability  to  produce  plasma  conditions 
approaching  those  of  the  natural  ionosphere.  Size  is  important  for  eliminating  wall  effects,  and 
parameters  approaching  ionospheric  LEO  values  are  important  for  relevance  to  space  concerns. 
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ChamberEnvironment 


The  SPSC  is  a  1.8-m  diameter  by  5-m  long  cylindrical  vacuum  chamber  shown  in  Figures  1  and 
2  which  can  be  evacuated  to  a  background  pressure  near  10'6  torr. 


Figure  1 


Coils  arranged  in  an  approximate  Helmholtz  configuration  provide 
an  axial  magnetic  field  variable  up  to  50  G.  There  is  the  capability 
to  pulse  the  field  to  kilogauss  levels.  Multiple  ports  and  an  internal 
3D  positioning  system  are  available  for  diagnostic  access.  The 

chamber  is  equipped  with  two 
microwave  plasma  sources 
[Walker  et  al.,  1994;  Bowles  et 
al.,  1996]  designed  with 
different  output  areas  for 
density  variation  and  differing 
experimental  requirements.  The 
wide  area  microwave  plasma 

source  is  shown  in  Figure  3.  In  addition  there  are  filament 
sources  which  can  provide  densities  as  high  as  1012  cm'3.  A 
variety  of  experiments  simulating  space  plasma  phenomena 
have  been  performed  in  the  chamber  [Walker  et  al.,  1995, 
Figure  3  i997;  Amatucci  et  al,  1996,  1998a, b] 

Table  I  outlines  the  basic  plasma  parameters  for  a  50-cm  diameter  argon  plasma  column. 
The  parameters  are  close  to  the  LEO  environment  so  that  scaling  concerns  in  this  chamber  are  a 
minimum.  While  the  ion-electron  temperature  ratio  is  lower  than  in  the  ionosphere,  the  base  ion 
temperature  does  correspond  well  to  that  in  the  ionosphere.  Also,  changes  in  operating 
conditions  for  the  microwave  plasma  source  (such  as  operation  under  higher  background 
pressures)  are  expected  to  bring  the  ratio  close  to  unity. 


Figure  2 
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Table  1.  SPSC  Plasma  Parameter  Regimes 


plasma  density 
neutral  density 
electron  temperature 
ion  temperature 
Debye  length 

electron  plasma  frequency 

electron  thermal  speed 

ion  thermal  speed 

ion  acoustic  speed 

electron  gyroradius 

electron  gyrofrequency 

ion  gyroradius 

ion  gyrofrequency 

ion-neutral  collision  frequency 

electron-neutral  collision  frequency 


106£  Nc<  1012  cm 3 
3  x  10"  £  N„  £  3  x  1014  cm'3 
Tc  -  0.5  eV 
T,  =  0.05  eV 
0.5  s  A.d<  .5  x  10 4  cm 
10  MHz£  fp,  £  300  MHz 
Vlc=  3  x  107  cm  s'1 
V,i=  5  x  104  cm  s'1 
Cs=l.l  x  105  cm  s'1 

pc  =  0. 1 6  cm  (B=10  G)  0.04  cm  (B=40  G) 

fcc=28  MHz  (  “  )  1 12  MHz  (  “  ) 

Pi  =  14  cm  (  “  )  3.6  cm  (  “  ) 

fcc=380  Hz  (  “  )  1.5  kHz  (  “  ) 

0.03 £  v,n/Q,  £30  (  “  )  0.01  £vin/Q,£8  (  “  ) 

4x104£vct/Qc£0.07(  “  )  104£Vc„/Qc£0.02(  “  ) 


Electronics:  charging  and  discharging  scenario 

A  typical  cycle  time  for  the  experiment  is  <  30  seconds.  This  includes  capacitor  charging, 
firing,  discharging  and  data  collection.  The  primary  analysis  routines  are  LAB  VIEW  software- 
based  and  AD  interfaces  employ  an  array  of  CAM  AC  digitizers. 


Figure  4 


The  SPEAR  III  program  included  testing  of  critical 
components  in  a  'mock-up'  (as  close  to  flight 
conditions  as  possible)  configuration  in  the  NASA- 
Lewis  Research  Center's,  Plum  Brook  Station  B2 
vacuum  facility.  The  mock-up  had  the  advantage 
of  allowing  a  direct  measurement  of  the  voltage 
difference  between  the  payload  body  and  the 
vacuum  chamber  wall.  This  was  carried  over  into 
the  NRL  experiment.  The  charging  electronics 
used  in  the  laboratory  testing  for  the  SPEAR  m 
program  were  brought  to  NRL,  and  the  initial 
configuration  was  designed  to  first  test  neutral  gas 
release  effects  from  a  single  cylinder,  as  shown 
schematically  in  Figure  4.  Also  shown  in  Figure  4 
are  the  relative  positions  of  the  light  sensitive  Figure  5 

diodes  used  in  an  investigation  of  the  dynamic 

sheath  expansion  during  discharge.  Details  on  the  charging  method  can  be  found  in  papers 
relating  to  both  the  SPEAR  I  [Allred  et  al,  1988;  Katz  et  al.,  1989]  and  SPEAR  HI  [Mandel  et 
al.,  1998;  Berg  et  al,  1995]  rocket  experiments.  For  the  mock-up  and  the  flight  experiments  one 
side  of  a  high-voltage  capacitor  on  the  payload  was  charged  to  +10  kV  and  the  other  side  was 
connected  through  a  relay  to  a  large  sphere  deployed  on  a  high-voltage  boom.  This  configuration 
was  chosen  so  that  in  a  plasma  with  the  characteristics  of  the  normal  F-region  ionosphere  the 
sphere/payload  combination  will  float  (with  respect  to  the  plasma)  with  at  least  -2  kV  appearing 
on  the  payload  with  the  remaining  voltage  appearing  on  the  sphere.  Thus,  as  the  voltage  on  the 
capacitor  decays,  the  payload  voltage  will  go  through  the  entire  range  of  voltages  normally 
associated  with  natural  charging,  i.e.,  0  to  -2  kV.  A  resistor  across  the  charging  capacitor  allows 
for  a  1  sec  decay  time  of  the  voltage.  Any  current  flowing  through  the  plasma  will  shorten  this 
decay  time.  During  the  experiments  the  decay  time  is  typically  between  0.3- 1.0  sec,  except  when 
a  large  breakdown  occurs  and  generates  a  high  current  discharge  path.  The  capacitor  is  then 
discharged  in  milliseconds.  For  the  rocket  experiment  the  charge  and  decay  cycle  was  set  to 
repeat  every  5  seconds.  The  SPSC  testing  was  done  on  a  single  'shot'  basis  where  the  charging 
and  data-taking  cycle  was  controlled  manually. 


For  the  laboratory  experiments  the  high  voltage  capacitor  could  be  connected  between  two 
cylinders,  or  a  single  cylinder  could  be  biased  with  respect  to  the  chamber  walls.  For  the 
experiment  described  here  we  chose  the  single  cylinder  configuration  and  applied  voltages  as 
high  as  -3  kV.  In  this  manner,  the  initial  SPSC  experimentation  eliminates  the  sphere  and 
applies  the  charging  voltage  directly  across  the  “payload”.  Discharging  of  the  simulated  payload 
in  the  SPSC  depends  on  the  current  collection  area  available  for  ions.  Shown  in  Figure  5  is  the 
initial  setup  of  the  two-cylinder  configuration  in  the  laboratory,  whereas  Figure  4  shows  the 
schematic  representation  of  the  initial  one  cylinder  experimental  arrangement.  In  the  laboratory 
simulation  the  small  cylinder  in  Figure  5  is  removed  from  the  experiment,  a  neutral  gas  release 
valve  is  placed  in  the  large  cylinder  and  kilovolt  level  charging  occurs  before  release.  The  single 
aluminum  cylinder,  which  is  10  cm  in  diameter  and  10  cm  in  length  and  was  mounted  in  the 


5 


center  of  the  chamber  is  shown  along  with  the  enclosed  puff  valve  system  schematically  in 
Figure  6a.  The  gas  supply  which  provides  the  neutral  gas  release  is  fed  to  the  plenum  inside  the 
cylinder  through  a  1/4  inch  plastic  line,  also  visible  in  Figure  6.  The  reservoir  is  connected  to  a 
solenoid  valve  so  that  when  the  valve  is  opened  the  gas  flows  through  and  into  the  throat  of  the 
supersonic  nozzle.  Two  separate  nozzles,  shown  in  Figures  6b  and  6c,  have  been  tested  so  far 
and  they  are  designed  to  produce  mach  numbers  of  Mach  3  and  Mach  9.  As  mach  number  is 
related  to  the  ratio  of  exit  area  to  throat  area  (shown  below),  the  separate  nozzles,  even  though 
machined  identically  at  the  throat,  have  different  lengths  from  throat  to  exit  and  thus  diverge  to 
different  exit  areas. 

During  the  discharges,  cylinder  voltage  and  current  collected  by  the  cylinder  were  measured 
as  functions  of  time.  A  typical  data  set  is  shown  in  Figure  7.  For  measurement  purposes,  several 
in-house-developed  diagnostics  were  used  including  fast-sweep  Langmuir  probes  [Siefring  et  al, 
1997b]  to  measure  electron  temperature,  emissive  probes  to  measure  plasma  potential,  photo 
diodes  to  monitor  neutral  excitation  and  double  probes  to  monitor  particle  flow  and  electric  field. 
In  addition,  a  high  voltage  sensor  [Siefring  et  al.,  1995]  was  developed  and  used  in  the 
laboratory  testing  in  the  determination  of  potential  at  various  positions  in  the  plasma. 


Figure  6a 


Figure  6b 


Figure  6c 
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Current  vs  Time 


Figure  7 


III.  Experimental  Results 

In  these  studies  we  have  done  work  related  to  the  effects  of  NGR  in  a  number  of 
discharges  where  the  voltage  applied  to  the  cylinder  is  varied  from  about  -600  V  to  -2400  V.  The 
gas  valve  is  typically  opened  a  short  time  (  « 10  ms)  after  charging  and  there  is  a  subsequent  local 
neutral  pressure  increase  by  several  orders  of  magnitude  in  the  near-space  of  the  release  valve. 
The  discharge  is  initiated  and  the  cylinder  voltage  falls  to  approximately  -200  V  as  the  collected 
current  rises  to  the  series  resistor-limited  value  of  about  1  amp.  This  behavior  is  seen  in  the  first 
two  plots  of  Figure  7  above. 
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minimum  potential  for  breakdown  (kV) 


Orrm  v  Cufranl 


It  should  be  noted  that  the  abrupt  rise  in  the 
cylinder  voltage  near  0.036  sec  is  an  artifact  and  not 
an  important  physical  effect.  Also,  much  of  the 
accompanying  experimental  noise  has  been  filtered 
from  this  data.  As  the  gas  breakdown  process 
continues,  electrons  are  accelerated  away  from  the 
cylinder  and  ions  are  collected.  Secondary  electrons 
are  created  and  also  accelerated  away  from  the 
cylinder  creating  further  neutral  excitation  and  other 
electron-ion  pairs.  The  light  emission  caused  by 
this  excitation  of  the  neutrals  was  monitored  by  an 
array  of  photo-diodes  arranged  in  increasing 
distances  from  the  discharge  region,  as  seen 
schematically  in  Figure  4.  The  placement  of  these 
detectors  is  the  first  stage  in  determining  expanding 
sheath  characteristics  during  the  discharge.  Typical 
traces  of  this  light  intensity  as  a  function  of  time 
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Figure  8 


Te  as  a  function  of  pulse  duration 


f  \ 


and  diode  number  are  shown  in  Figure  8  with  the  intensity  seen  to  fall  off  as  a  function  of 
distance.  The  use  of  band  pass  filters  over  the  photodiode  demonstrates  that  the  photons  result 
primarily  from  neutral  collisions. 

In  Figure  9  the  electron  temperature 
measurements  as  a  function  of  pulse  duration  on  each  of 
three  probes  at  varying  axial  positions  from  the  discharge  3 
region  show  enhanced  electron  temperature  at  the 
farthest  probe  position.  This  enhancement  is  consistent 

with  the  falloff  in  2 

min.  breakdown  potential  vs.  plenum  pressure  nCUtral  density  and  the  2 

inelastic  cooling  rate  | 
of  the  electrons.  The 
discharge  duration  is  1 

typically  bounded 
initially  by  the  release 
of  neutral  gas  and 
finally  by  the  decrease 
of  potential  below  the 
minimum  necessary  to 

sustain  ionization.  Q 

The  gas  puff  duration  igure 

is  significantly  longer  than  the  discharge  time. 

Figure  10  is  a  plot  of  the  minimum  potential  necessary 
Figure  10  t0  *n't*ate  breakdown  versus  plenum  pressure  for  three  separate 
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gas  puff  species  at  a  single  fixed  mach  9.  Each  point  in  Figure  10  corresponds  to  an  average  of 
10  separate  experimental  runs.  From  this  plot  one  can  conclude  that  krypton  and  argon  initiate 
breakdown  at  a  lower  potential  than  neon  for  the  same  gas  plenum  pressure,  suggesting  that 
perhaps  heavier  neutral  gases  initiate  breakdown  at  lower  potentials  (See,  however,  the  sections 
on  theory  and  comparison  to  experiment  below).  Recent  models  [Langely  et  al.,  1984;  Dietz 
and  Sheffield,  1975]  allude  to  a  velocity  dependence  on  secondary  electron  emission  that  would 
seem  to  favor  lighter  gases.  In  Section  IV  we  present  a  breakdown  theory  which  addresses  these 
possibilities  and  describes  the  data  fairly  well.  We  are  in  the  process  of  building  on  the  initial 
study  by  investigating  the  effectiveness  of  other  gas  species  for  mitigation  of  spacecraft  charging. 

With  a  goal  toward  minimizing  gas  usage,  and  therefore  storage,  we  have  studied  the 
breakdown  characteristics  versus  puff  gas  flow  rate.  Figures  11  and  12  show  photographs  of  the 
discharges  in  argon  and  neon  for  the  mach  3  and  mach  9  nozzles  separately.  In  general,  the 
discharge  appears  more  elongated  for  the  higher  mach  number. 


Argon  puff 


Neon  puff 


Figure  11 


Figure  12 


min.  breakdown  potential  vs.  gas  mass  flow  rate 
argon,  B  =  1  O  G,  background  plasma  density  n=  1  Ob  cm  3 


Figure  13  shows  initial  results  of 
the  discharge  studies  for  Argon  only. 
Each  data  point  is  again  an  average  of 
10  separate  shots.  The  preliminary 
conclusion  from  this  data  is  that  the 
mach  9  nozzle  valve  discharges  the 
cylinder  at  lower  potentials  than  the 
mach  3  valve.  We  are  preparing  further 
experimentation  in  this  area  with  a 
wider  combination  of  release  gases  and 
valves. 


0.0  0.2  0.4  0.6  O.e  1.0 


gas  mass  flow  rate  (g/sec) 

Figure  13 
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IV.  Theory 


The  Discharge  Model 


The  study  of  neutral  gas  breakdown  in  the  presence  of  high  potentials  has  been  underway 
since  the  beginning  of  this  century  [Townsend,  1900,  1901;  Loeb,  1939;  Raizer,  1997],  The 
fundamental  assumption  is  that  the  number  of  electrons  grows  exponentially  in  an  avalanche 
fashion  as  one  proceeds  from  cathode  toward  anode: 


X 

N(x)  W()exp[  J  a(x)dx] 
o 


(1) 


where  N0  is  the  number  of  electrons  released  from  the  cathode  and  a  is  the  Townsend  ionization 
coefficient  and  is  defined  as  the  number  of  ionization  events  performed  along  a  1  cm  path  in  the 
electric  field  direction.  Equivalently, 


(2) 


with  Vj  the  ionization  frequency  and  vd  the  electron  drift  speed,  a  is  therefore  a  function  of  the 
neutral  density  (or  pressure)  and  the  electric  field.  An  empirical  formula  for  a  , 


a2n 

a=AlNe  E 


(3) 


is  used  in  much  of  the  theoretical  and  numerical  analysis.  Here  A,  and  A2  are  determined  from 
the  experimental  data  for  a  particular  gas  and  N  is  the  neutral  gas  density.  (We  note  here  that 
often  this  expression  is  given  in  terms  of  neutral  gas  pressure.  This  causes  complications, 
however,  when  temperature  varies  and  the  usual  standard  room  temperature  is  not  the  applicable 
temperature.  To  use  pressure  in  the  conventional  sense  one  has  to  make  the  distinction  that  the 
pressure  given  is  related  to  the  neutral  density  through  standard  room  temperature,  and,  in  any 
case,  this  is  confusing.  For  this  reason  we  prefer  to  express  the  relation  through  N.)  For  inert 
gases,  however,  a  better  fit  to  the  experimental  data  is  provided  by  Raizer,  [1997], 
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(4) 


-b2(|)2 

a=BlNe  E 


where  Bj  and  B2  are  empirical  “constants”  determined  from  the  ionization  curves. 

The  condition  for  a  self-sustaining  discharge  [Raizer,  1997]  is, 

OO 

I=Jcbc  a(x)  z  ln(l+Yi*‘)  (5) 

o 


where  Yi  is  the  coefficient  for  electron  emission  from  the  cathode  per  incident  ion.  For  Yi  ~ 

0.15,  for  example,  the  right  hand  side  of  this  equation  equals  2.  The  minimum  voltage  should  be 
100  V  or  more,  based  on  the  data  for  Townsend  discharges  (space-charge  free)  and  glow 
discharges  (with  space  charge).  At  this  voltage,  the  gas  density,  N,  times  the  gap  separation,  d, 
ranges  from  7  x  101S  to  5  x  1016  cm'2,  depending  on  the  gas,  the  cathode  material  and  the  spatial 
profile  of  E/N.  At  higher  or  lower  values  of  Nd,  larger  voltages  are  needed  to  produce 
breakdown.  This  suggests  that  a  gas  puff  can  discharge  an  object  of  characteristic  radius  R, 
provided  it  is  charged  (negatively)  to  a  few  hundred  volts  or  more  and  NR  lies  between  7  x  10 15 
and  5  x  1016  cm"2.  Here  N  is  the  gas  density  averaged  over  R.  In  the  present  experiments,  R~  10 
cm,  which  suggests  that  the  optimal  value  of  N  lies  between  1015  and  1016  cm'3.  Such  densities 
are  achievable,  according  to  the  analysis  below.  Note  that  the  density  required  varies  inversely 
with  the  size,  R,  of  the  object. 

Analysis  of  Eqs  (4)  and  (5)  shows  that  the  minimum  voltage  for  breakdown,  <j)min,  occurs 
near  the  peak  of  the  ionization  curve,  where  a/N  =  (a/N)max  and  E/N  =  (E/N)min.  Setting  (E/N)min 
=  (J)min/NR  and  amaxR  equal  to  In  (1+Yi'1)  based  on  condition  (2),  we  obtain 

<J>mi„=  [(E/N)min/(a/N)max]ln  (1+Yi  -1)-  In  noble  gases  (a/N)max  ~  3  xlO'16  cm2  and  (E/N)min  ~  3  x 
10'14  V-cm2.  The  minimum  breakdown  voltage  is  then  <J)min  ~  100  In  (1+Yi '')  ~  200  V.  Note  that 
<j)min  scales  directly  with  the  right  hand  side  of  Eq.  (5).  It  might  be  noted  that  this  voltage  is 
roughly  an  order  of  magnitude  smaller  than  the  minimum  voltage  for  self-sustained  streamers  as 
given  by  Raether,  [1939] 


Neutral  Gas  Release  (NGR)  -  Description  of  the  Expansion 


Equations  (3)  and  (5)  indicate  that  two  parameters  determine  whether  breakdown  occurs 
for  a  particular  gas.  The  first  is  the  field  profile,  E(z),  which  is  determined  by  the  charging 
voltage  and  system  geometry.  The  other  is  the  gas-density  profile,  N(z)  =  F(/v(z)A(z),  where  F0 
is  the  total  flux,  v(z)  is  the  local  flow  velocity,  and  A(z)  is  the  local  transverse  area  of  the  gas 
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puff.  We  compute  F0,  v(z)  and  A(z),  relying  in  part  on  the  early  work  of  Barrere,  [I960], 
Consistent  with  this  approach,  we  assume  isentropic  expansion  and  isotropic  temperature  along 
z.  In  addition  we  point  out  some  applicable  features  associated  with  definition  of  nozzle  mach 
number. 


(i)  No  nozzle:  Free  expansion  of  an  ideal  gas 


Consider  a  large  pressurized  box  with  a  small  hole  in  one  side.  Particles  inside  are 
assumed  to  have  a  Maxwellian  velocity  distribution,  and  only  those  with  vz  >  0  can  escape 
through  the  hole.  If  the  expansion  is  into  complete  vacuum,  then  the  mean  forward  velocity  is 
given  by. 


/' 


dvve 


2  r. 


<v,> 


IT 

c 

\  nm 


f 


dve 


2  T, 


(6) 


where  Tc  is  the  chamber  temperature  and  m  is  the  molecular  mass.  The  mean  forward  speed  is 
thus  less  than  0.7vs,  where  vs  =  (yTc/m)1/2  is  the  sound  speed  and  y  £  1.2  is  the  adiabatic 
constant.  In  passing  through  the  hole  the  transverse  temperature  is  unchanged,  Tx  =  Tc,  while  the 
longitudinal  temperature  falls  to  Tz  =  <mvz2>  -  m<vz>2=(l-2/7i)Tc  ~  0.36  Tc.  Gas  collisions  will 
eventually  restore  temperature  isotropy  and  in  a  noble  gas  Tz  and  T±  approach  0.8  Tc.  Thus  only 
twenty  percent  or  so  of  the  initial  thermal  energy  is  realized  as  forward  energy. 


(ii)  Expansion  with  nozzles 


However  a  much  higher  conversion  efficiency  is  obtainable  using  a  gas  nozzle.  In  a 
nozzle  the  flow  is  primarily  one  dimensional  through  an  area  which  varies  with  distance.  The 
area  first  constricts  to  a  minimum  at  the  throat  where  the  flow  speed  becomes  the  local  sound 
speed,  v,=(yTt/m)l/2  where  T,=(2/y+l)Tc  is  the  gas  temperature  at  the  throat.  Past  the  throat  the 
nozzle  area  again  increases  as  can  be  seen  in  the  nozzle  diagrams  of  Figure  6.  Using  part  of  the 
analysis  by  Barrere  [1960,  Ch.  2  ],  we  can  show  that  the  nozzle  area  A  and  velocity  v  are  related 
to  the  values  at  the  throat  by, 
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(7) 


£ 

zf 


v  r'Y  +  l  vi  _/Y_1 


-  -[(^XI-C^X-)2)] 

v  2  Y+1  v, 


i 

Y-l 


•This  equation  indicates  that  v/vt  is  double  valued  for  a  given 
value  of  A.  This  relation  is  plotted  in  Figure  14  which  shows 
that  in  the  constricting  region  before  the  throat ,  v<vt;  and  in 
the  expanding  region  after  the  throat,  v  >  vt.  The  velocity 
thus  rises  as  the  nozzle  constricts  and  as  it  expands.  The 
velocity  at  the  nozzle  exit  depends  on  the  final  area  ratio  but 
is  always  limited  (for  A  -  °° )  to 


Figure  14 


v„  <  v 

e  max 


y  +  1 

N  Y-l 


v'  =  N 


Y-l 


(8) 


Nozzle  mach  numbers  are  typically  quoted  by  using  the  ratio  of  ve  to  the  local  sound 
speed,  c  s„  at  nozzle  exit  or, 


cs\ 


Jy-L 


(  l)2(Y-0[  e  ]2(Y+D 

c,  r  'a. 


(9) 


Here  Ae  is  the  exit  area  and  Te  is  the  exit  temperature  which  we  can  find  from, 

T  2  A.  v. 

T  1+y  Av/ 

c  *  e  e 


(10) 


where  the  coefficient  T'  is  given  by, 


r' 


Y*! 

2(y-1) 


(11) 


13 


The  flow  rate  F0  equals  the  mass  flow  rate  divided  by  the  molecular  mass,  m. 
Alternatively,  F0  can  be  computed  from  the  chamber  pressure  pc  and  the  throat  area  A,  using, 


F0  =  where  T  =  yfiV 


(12) 


(iii)Expansion  after  leaving  nozzle 


We  next  determine  the  velocity  v(z)  and  the  area  A(z)  outside  the  nozzle.  Because  there 
is  residual  energy  in  the  gas,  the  expansion  accelerates  longitudinally  and  transversely.  The 
longitudinal  acceleration  is  negligible,  however,  because  the  exit  temperature  is  small  compared 
with  the  chamber  temperature,  Te«Tc.  For  example,  for  the  M3  nozzle,  Te/Tc  =  0.05  according 
to  Eq.  (10);  for  the  Mach  9  nozzle  this  ratio  is  even  smaller.  Therefore,  past  the  nozzle  it  is 
legitimate  to  approximate  v(z)  ~  ve. 


Transverse  acceleration  is  generally  small  as  well.  To  show  this,  we  observe  that  the  gas 
can  be  treated  as  a  jet  moving  with  a  fixed  longitudinal  speed,  ve,  outside  the  nozzle.  Since  the 
transverse  velocity  is  presumed  small,  vA«ve,  the  jet  can  be  treated  as  a  freely  expanding,  long 
beam.  In  that  case  the  area  expands  as  [  Fernsler  et  al.,  1994] 

IT 

A(z)  =  Ae+  2zJtzA~  tan0o  +  7tz2[tan20o+ — e-\  (13) 

mv~ 

e 


where  0O  =  v/ve  at  z=0  is  the  half-angle  of  the  nozzle. .  This  expansion  produces  a  larger  area  at 
a  given  z  than  simple  cylindrical  expansions  due  to  the  presence  of  thermal  energy. 


V.  Comparisons  of  Theory  and  Experiment 


Breakdown  potentials  versus  gas  flow 


Using  the  criterion  outlined  in  Eq  (5)  above  we  plot  in  Figures  15  through  17  the 
predicted  breakdown  voltages  for  the  three  gases  used  in  the  experiment  and  for  each  of  the  two 
nozzles.  The  plot  of  Figure  1 5  also  shows  the  argon  data  superimposed  on  the  theoretical  plots. 
From  these  plots,  the  lightest  gas  tested,  neon  (Figure  17),  requires  significantly  higher  levels  of 
release  gas  to  break  down  for  a  given  cylinder  charging  voltage  than  either  argon  or  krypton. 
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This  result  is  consistent  with  the  data  shown  in  Figure  10  where  breakdown  potential  versus 
plenum  pressure  is  plotted  for  each  of  the  three  release  gases  using  the  mach  9  nozzle. 
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If  we  examine  the  breakdown  condition  of  Eq.(5)  in  somewhat  more  detail  the  integral  appears 
as, 


/  - 


=  f  N{z,p)c] 


/ 


E(z) 


e  *  N{:  p)dz 


(14) 


where  N(z,p)  is  neutral  density  given  by, 


N(z,p) 


Fo(P) 

Mz)ve  ’ 


(15) 


and  zL  is  the  lower  integration  limit.  F0(p)  is  the  conserved  neutral  particle  flux  from  the 
plenum,  A(z)  is  as  given  in  Eq.  (13)  and  ve  is  the  exit  velocity  which  is  calculated  from  the  area 
ratio  of  Eq.  (7).  This  relation  assumes  that  the  velocity  of  the  exit  gas  is  constant  as  a  function  of 
distance  from  the  release  point.  We  note  specifically  that  zL  corresponds  to  the  surface  of  the 
charged  object  where  spherical  symmetry  is  assumed,  i.e.,  z  =  0  is  taken  as  the  origin  of  a 
charged  sphere  for  definition  of  the  spatial  falloff  of  the  potential.  In  this  analysis  then,  the 
electric  field  and  potential  have  the  simple  one-dimensional  form, 

w)  -  -yj-)  ••  -  yjr\)  (i6) 

2  Z2 


where  Vch  is  the  charging  potential.  With  the  definition  of  zL  as  above,  we  are  assuming  that  the 
exit  area  of  the  nozzle  is  at  the  surface  of  the  cylinder;  therefore  A(zL)  corresponds  to  Ae  in  Eq. 
(7)  and  N(zL,p)  is  the  gas  density  at  the  exit  area  face  (note,  not  the  throat).  F0(p)  is  dependent 
upon  the  pressure,  p,  and  temperature,  T,  in  the  plenum,  the  release  gas  molecular  mass,  p,  and 
the  area  of  the  nozzle  throat,  A, ,  i.e., 


Fo  (P) 


r 


p  At 

sf^Rf 


(17) 


In  this  equation,  the  molecular  mass,  m,  is  in  kg,  p  in  pascals,  A,  in  m2,  T  in  Kelvin,  T  is  defined 
in  Eq.  (12),  and  R  is  the  universal  gas  constant. 
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Breakdown  dependence  on  release  gas  mass 


Since  the  density,  as  seen  in  Eq.  (15),  appears  in  the  breakdown  integral  Eq.  (14)  both  as 
a  multiplicative  factor  and  also  in  the  exponent  under  a  square  root,  it  is  not  immediately  obvious 
what  overall  effect  an  increase  of  mass  produces  without  numerical  evaluation  of  the  integral  as 
we  have  done  here.  In  fact,  however,  we  can  show  that  Ncm(z,p)  is  independent  of  the  release  gas 
mass  and  therefore  that  the  breakdown  integral  of  Eq.  (15)  is  largely  determined  by  the  ionization 
coefficients  for  the  various  gases;  therefore,  although  the  breakdown  is  dependent  upon  the  gas 
species  through  these  coefficients,  it  is  not  a  function  of  the  mass.  In  fact  we  can  show, 

m  =  4hNe  (18) 


where  Ne  is  the  gas  density  at  nozzle  exit  and  is  related  to  density  in  the  chamber  through, 


T  — 

N,  -  ( fy-'Nc  (19) 


Eq  (18)  results  from  simple  flux  conservation  and  the  assumption  of  (nearly)  constant  flow 
velocity  outside  the  nozzle.  Eq.  (19)  arises  from  Eqs.  (6)  and  (7)  and  the  relation  between 
velocity  in  the  nozzle  throat  and  sound  speed  in  the  plenum, 


v 


t 


(20) 


The  independence  of  density  and  therefore  of  breakdown  on  gas  mass  is  apparently  the 
explanation  of  the  experimental  result  shown  in  Figure  10  that  the  breakdown  potential  gas  flow 
requirement  for  krypton  and  argon  is  similar  whereas  the  difference  between  these  two  and  neon 
is  substantial.  This  occurs  even  though  there  is  approximately  a  factor  of  two  mass  difference 
between  each  mass  and  so,  at  first  glance,  one  might  be  inclined  to  guess  a  more  or  less  linear 
dependence  on  release  gas  mass.  As  covered  above,  this  dependence  is  largely  a  function  of  the 
ionization  curves  for  the  different  species  and  not  the  mass. 


Breakdown  dependence  on  release-gas  mach  number 
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As  demonstrated  above  the  gas  density  at  a  fixed  distance  from  the  nozzle  is  independent 
of  mass;  however,  not  of  exit  velocity  and  therefore  (local)  mach  number  as  seen  in  Eq.(15). 
Since  the  variation  of  pressure  in  the  plenum  is  essentially  an  isothermal  process  there  is  no 
effect  on  the  gas  exit  velocities  from  this  variation.  There  is  only  an  adjustment  in  gas  density 
which  occurs  with  varying  plenum  pressure.  At  a  fixed  pressure  and  fixed  distance  from  the 
nozzle  therefore,  we  plot  the  density  dependence  on  the  ratio  of  exit  velocities  in  Figure  18.  We 

see  that  as  the  exit  area  increases  (or  as  the 
local  mach  number  or  velocity  ratio 
increases)  the  density  decreases.  Although 
not  shown  the  density  continues  to  decrease 
to  the  point  that  at  the  limit  that  the  velocity 
ratio  is  2,  we  produce  an  infinite  area  and 
hence  a  density  decrease  to  zero.  Figure  19 
shows  a  plot  of  the  breakdown  integral  as  a 
function  of  the  velocity  ratio  for  neon  at  a 
fixed  position  and  pressure.  The  value  of 
the  breakdown  integral  is  seen  to  be  a 
monotonic  function  of  the  velocity  ratio, 
and  hence  of  the  valve  mach  number,  up 
until,  the  limit  of  infinite  area  where  the 
function  decreases.  This  occurs  because  of 
the  steep  decline  in  the  gas  density  levels 
v  m  beyond  a  certain  distance  from  the  nozzle. 

Figure  18  These  theoretical  arguments  would  then 

suggest  that  the  higher  mach  number 

nozzles  should  produce  breakdown  with  less  gas  pressure  than  lower  mach  number  nozzles. 

These  results  are 
corroborated  by  the 
experimental  plots  as 
pointed  out  above. 
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Figure  19 
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VI.  Summary 


We  have  conducted  preliminary  laboratory  simulations  and  developed  a  basic  theoretical 
model  of  the  effects  of  neutral  gas  release  on  space  vehicle  charging  in  the  LEO  space 
environment.  The  experiments  were  performed  in  the  NRL  Space  Physics  Simulation  Chamber. 

The  results  presented  here  are  the  first  stage  of  a  more  exhaustive  investigation  into  the 
optimum  choice  of  a  release  gas  and  release  speed  in  order  that  spacecraft  which  acquire 
potentials  either  through  natural  or  artificial  means  can  discharge  these  voltages  at  the  lowest 
level  possible.  The  Townsend  discharge  model  and  breakdown  criterion  expressed  through  the 
integral  of  Eq  (5)  explain  and  predict  the  experimental  results  quite  well.  Because  of  the  success 
of  this  model  we  are  in  the  process  of  continuing  work  with  it  in  order  to  guide  certain  phases  of 
the  experimental  effort  and  to  suggest  nozzle  and  gas  release  designs  for  optimum  performance. 

At  the  present  stage  of  the  experiment  we  are  able  to  conclude  that  argon  and  krypton 
tend  to  discharge  the  payload  at  lower  potential  than  neon  and  that  the  higher  mach  number 
(M9)  valve  also  is  better  at  discharging  for  a  given  plenum  pressure  than  the  lower  mach  number 
valve  (M3).  We  are  continuing  this  experiment  with  studies  using  different  neutral  gases  and 
mixtures  of  neutral  gases,  with  different  cylinder  surfaces  and  with  varying  electron  and  neutral 
gas  environments.  In  addition  we  are  investigating  the  effect  of  the  orientation  of  the  puff  with 
respect  to  the  local  magnetic  field. 


Acknowledgment: 

We  gratefully  acknowledge  the  support  of  the  Office  of  Naval  Research.  In  addition,  the  Ballistic 
Missile  Defense  Organization  (BMDO)  was  responsible  for  early  support  in  the  SPEAR  effort 
from  which  much  of  the  pulsing  electronics  was  derived. 


19 


VII.  Figure  Captions 

Figure  1.  A  schematic  diagram  of  the  NRL  Space  Physics  Simulation  Chamber  showing 
diagnostic  equipment 

Figure  2.  A  photograph  of  the  NRL  SPSC  facility. 

Figure  3.  A  photograph  of  the  wide-area  microwave  plasma  source  showing  magnetron  and 
coupling  wave  guides. 

Figure  4.  A  schematic  diagram  of  the  charging  electronics  and  cylinder/nozzle  inside  the  SPSC 

Figure  5.  A  photograph  of  the  the  two-cylinder  charging  arrangement  inside  the  SPSC 

Figure  6.  Engineering  drawings  of  a.  the  charged  cylinder  with  cut-away  view  of  gas  plenum  and 
nozzle,  b.  the  mach  3  nozzle,  c.  the  mach  9  nozzle. 

Figure  7.  Plots  of  experimental  collected  current,  cylinder  voltage  and  impedance  as  a  function 
of  time  during  a  discharge. 

Figure  8.  Experimental  traces  of  current  collection  and  light  intensity  as  measured  by  diodes  at 
varying  distances  from  the  cylinder  during  discharge 

Figure  9.  Experimental  electron  temperature  as  measured  by  Langmuir  probe  sweeps  at  varying 
distances  from  the  cylinder  as  a  function  of  pulse  duration. 

Figure  10.  Experimental  plots  of  minimum  potential  necessary  for  breakdown  for  each  of  three 
gases  as  a  function  of  plenum  pressure.  Mach  9  (M9)  nozzle. 

Figure  11.  A  photograph  of  an  argon  gas  puff  (M3  nozzle)during  discharge. 

Figure  12.  A  photograph  of  a  neon  gas  puff  (M9  nozzle)  during  discharge 

Figure  13.  Experimental  plots  of  minimum  potential  necessary  for  breakdown  for  each  of  the 
two  nozzles,  M3  and  M9  as  a  function  of  gas  mass  flow  rate. 

Figure  14.  A  theoretical  plot  of  the  ratio  of  nozzle  throat  area  to  exit  area  versus  ratio  of  exit 
velocity  to  sound  speed  in  nozzle  throat.  This  theoretical  relation  is  plotted  from  Eq.  (7). 

Figure  15.  A  theoretical  plot  of  minimum  breakdown  voltage  versus  plenum  pressure  for  argon 
with  two  different  nozzles,  M3  and  M9.  Experimental  breakdown  for  the  M9  nozzle  is  plotted  as 
triangles  for  comparison.  As  derived  from  the  breakdown  integral  criterion  of  Eq  (5). 

Figure  16.  A  theoretical  plot  of  breakdown  voltage  versus  plenum  pressure  for  krypton  with  two 
different  nozzles,  M3  and  M9.As  derived  from  the  breakdown  integral  criterion  of  Eq  (5). 

Figure  17.  A  theoretical  plot  of  breakdown  voltage  versus  plenum  pressure  for  neon  with  two 
different  nozzles,  M3  and  M9.  As  derived  from  the  breakdown  integral  criterion  of  Eq  (5). 

Figure  18.  A  theoretical  plot  of  neon  neutral  gas  density  (Eq.  (15))  at  fixed  distance  from  the 
cylinder  and  constant  pressure  versus  release  gas  mach  number. 

Figure  19.  A  theoretical  plot  of  the  results  of  the  the  breakdown  integral  criterion  of  Eq.  (5)  at 
lixcd  distance  from  the  cylinder  and  constant  pressure  versus  release  gas  mach  number 
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